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J ' Abstract. We study the evolution of a weakly convex surface Sq in M'^ with 

C^ flat sides by the Harmonic Mean Curvature flow. We establish the short time 

existence as well as the optimal regularity of the surface and we show that 
the boundaries of the flat sides evolve by the curve shortening flow. It follows 

^^ ' from our results that a weakly convex surface with fiat sides of class C*^'^, for 

*vj ' some A: G N and < 7 < 1, remains in the same class under the flow. This 

distinguishes this flow from other, previously studied, degenerate parabolic 
equations, including the porous medium equation and the Gauss curvature 

P^ ' flow with flat sides, where the regularity of the solution for t > does not 

.^^ ' depend on the regularity of the initial data. 
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1. Introduction 



We consider the motion of a compact, weakly convex two-dimensional surface 
So in space K"^ under the harmonic mean curvature flow (HMCF) 

, ^ ! , , dP K 

\0 , (HMCF) = N 



dt H 

0^ ' where each point P of So moves in the inward normal direction N with velocity 

C^ I equal to the harmonic mean curvature of the surface, namely the harmonic mean 

^' ■ K \^ A2 



H Ai + A2 

(^ , of the two principal curvatures Ai, A2 of the surface. 

The existence of solutions to the HMCF with strictly convex smooth initial data 
was first shown by Andrews in [2]. He also showed that, under the HMCF, strictly 

^^ ' convex, smooth surfaces converge to round points in finite time. In [ID] , Dieter 

established the short time existence of solutions to the HMCF with weakly convex 
smooth initial data and mean curvature H > 0. More precisely. Dieter showed that 
if at time t = the surface So satisfies K > and H > 0, then there exists a 
unique strictly convex smooth solution Sj of the HMCF defined on < i < r, for 
some T > 0. By the results of Andrews, this solution exists up to the time where 
its enclosed volume becomes zero. However, the highly degenerate case where the 
initial data is weakly convex and both K and H vanish in a region is not studied 
in [10] . 

We will consider in this work the evolution of a surface Sg with flat sides by 
the HMCF. The parabolic equation describing the motion of the surface becomes 
degenerate at points where both curvatures K and H become zero. Our main 
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objective is to study the solvability and optimal regularity of the evolving surface 
for i > 0, by viewing the flow as a free-boundary problem. It will be shown that 
a surface Eo of class C'''^ with fc e N and 0<7<latt = 0, will remain in the 
same class for t > 0. In addition, we will show that the strictly convex parts of the 
surface become instantly C°° smooth up to the flat sides and the boundaries of the 
flat sides evolve by the curve shortening flow. 

For simplicity we will assume that the surface Sq has only one flat side, namely 
S = Si U S2, with Si flat and S2 strictly convex (both principal curvatures are 
strictly positive). We may also assume that Si lies on the z — plane and that S2 
lies above this plane since the equation is invariant under rotation and translation. 
Therefore, the lower part of the surface So can be written as the graph of a function 

z = h{x,y) 

over a compact domain i? C M^ containing the initial flat side Si. Let F denote 
the boundary of the flat side Si. We deflne g = h^ , for some < p < 1. Our main 
assumption on the initial surface Sq is that it satisfies the following non- degeneracy 
condition (*): 

(*) \Dg{P)\ > A and grriP) > \ for aU P G F 

for some number A > 0. Here t denotes the tangential direction to the level sets of 
g and grr denotes the second order derivative in this direction. 
Under the above conditions, our main results show that for t G (0,T): 

(1) The HMCF admits a solution St = (Si)t U (Sa)* of class C'''^ , for some 
fc G N and < 7 < 1 depending on p, which is smooth up to F^ = 9(Si)t. 

(2) (Si)t is flat and its boundary F^ evolves by the curve shortening flow. 

The fact that the solution Sf remains in the class C'^''^ distinguishes this flow from 
other, previously studied, degenerate free-boundary problems (such as the Gauss 
curvature flow with flat sides, the porous medium equation and the evolution p- 
laplacian equation) in which the regularity of the solution for t > does not depend 
on the regularity of the initial data. 

We define 6 to be the class of weakly convex compact surfaces Sg in R^ so that 
S = Si U S2, where Si is a surface contained in the plane z — and S2 is a strictly 
convex and smooth surface contained in the half-space z > 0. The main result 
states as follows: 

Main Theorem. Assume that at time i = 0, Sq is a weakly convex compact 
surface in M^ which belongs to the class & so that the function g — h^ defined as 
above is smooth up to the interface F and satisfies fk). Then, there exists a time 
T > such that the HMCF admits a solution St G 6 on [0,T). Moreover, the 
function g ^ h^ , defined as above for Sj, is smooth up to the interface z — and 
satisfies (-k) on (0, T] . In particular, the interface Ft between the flat side and the 
strictly convex side is a smooth curve for all t in < t < T and it evolves by the 
curve shortening flow. 

Sketch of the proof. A standard computation shows that when Sj solves the HMCF, 
the function h evolves by the equation 

h h — /)2 

^'' '^ {l + hl)h,,~2h,hyh,y-il + hl)hyy "^ ^ ^ ' 
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The HMCF can be seen as a free boundary problem arising from the degeneracy 
near the fiat side of the fully nonlinear parabolic PDE which describes the flow. We 
will show in section [5] that, via a global change of coordinates, this free boundary 
problem is equivalent to an initial value problem on Px [0, T], with D = {{u, v); u^ + 
v'^ < 1}, namely 

, , f Mw^O onVx [0,T] 

\ w — wq at t ^ 

The operator Mw = wt — F{t,u,v,w, Dw, D'^w), is a fully non-linear operator 
which becomes degenerate at dV, the boundary of V. We will apply the Inverse 
Function Theorem between appropriately defined Banach spaces to show that this 
problem admits a solution. 

The linearization of the operator M at a point w close to the initial data wq, 
can be modeled, after straightening the boundary, on the degenerate equation 

(1-3) ft = Z^ ail/^^ + 2 Z ai2fzy + 022 fyy + 6l Z /z + ^2 fy 

on z > with no extra conditions on / along the boundary z = 0. We observe 
that the diffusion in the above equation is governed by the Riemannian metric 

ds^ — ds^ -\~\dt\ where 

dz^ 
df = -^+dy'^. 

z^ 

The distance (with respect to the singular metric s) of an interior point (z > 0) 
from the boundary (z = 0) is hence infinite. This fact distinguishes our prob- 
lem from other, previously studied, degenerate free-boundary problems such as the 
degenerate Gauss curvature flow and the porous medium equation. 

The plan of the paper is the following: in section [5] we will introduce a local 
change of coordinates that fixes the free-boundary T in equation (|l.ip . We will 
compute the linearization of equation (|l.ip in this new change of coordinates, and 
motivate the use of the appropriate Banach spaces C^+^-p for our problem. The 
detailed definition of these Banach spaces will be given in section [31 where we will 
also present the appropriate Schauder estimates for our problem. In section S] we 
will study the fully-nonlinear degenerate equations (|1.2p and establish the short 
time existence for such equations in the Banach spaces C^ +"'''. The global change 
of coordinates and the proof of the Main Theorem will be given in sections [5] and 
ini respectively. In the last section we will establish the comparison principle for 
equation (jl.ip and characterize our solutions as viscosity solutions. 

Acknowledgments. The authors wish to thank G. Huisken for suggesting the 
problem and R. Hamilton for many stimulating discussions. 

2. Local Change of Coordinates 

We will assume throughout this section that the surface Eq belongs to the class 
6. Let St be a solution to the HMCF on [0,T), for some T > in the sense that 
St = (Si)t U (^2)4, with (Si)t flat and (S2)t strictly convex. Let Pq (ccq , yo , ^o ) be 
a point on the interface Ttg, for to > sufficiently small. Then, the strictly convex 
part of surface (^2)45, t < to can be expressed locally around Pq as the graph of a 
function z = h{x, y, t). Let g be defined hy g = hP, for < p < 1, such that : g is 
smooth up to the interface and satisfies condition (•). We can assume, by rotating 
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the coordinates, that at the point Pq the normal vector to Ff,, facing outwards the 
flat side S(io) is parallel to the x-axis, so that at Pq we have 

g^{Po) > and gy{Pn) = 0. 

Then we solve locally around the point Pq the equation z ~ h(x, y, t) with respect 
to X. This yields to a map x — f{z, y, t). The condition on g expressed in terms of 
/ gives the following non- degeneracy condition (irk) : 

\ ^ Jzy ~ Jyy J 

in the sense that both eigenvalues of the above matrix are bounded from below by 
a number A > 0. 

Since / is the inverse of h and the HMCF is invariant under rotation, the function 
/ satisfies the same equation as ft, on z > ; 

^' ' ^* " (1 + /2)/.. - 2 f^fyf^y - (1 + f^)fyy ^ 

We will construct a smooth solution to this equation by using the Inverse Function 
Theorem. To do so, we will define the Banach space C'^+°''P in the next section. 
According to our notation, the constants a and p indicate "how the surface becomes 
flat" , while s refers to the hyperbolic metric which governs the problem. 
We will prove in the next sections that when / € C^+"'P and satisfies condition (**), 
then the equation (j2.ip becomes degenerate at z = implying that: 



Jyy 



(2.2) ft = ^^ at the interface z = 0; 



This is equivalent to say that the free boundary Ft evolves by the curve shortening 
flow. 

3. The c'^+°''P space and Schauder estimates 

Let ^ be a compact subset of the half space { (z, y) G M^ : z > } such that 
(0, 0) e A. Then, we define: 

A° := {y e M :(0,y) e A} 

A := {{w,y) eM.'^ ■.w = \nz, {z,y) e A, zj^ 0} 

Qt := Ax [0,r], T>0 

Q°j. := A° X [0, T] 

Qt ■■= Ax[0,T]. 

Let < p < I. Given a function / on ^ wc define: 

riy) ■■= /(0,y) 

f{w,y) := e-P^{f{z,y)^r{y)) 

with w = Inz, for z > 0. 

Analogously, given a function / on Qt we define: 

ny,t) := f{0,y,t) 

fiw,y,t) := e-P-(/(z,y,t)-r(y,i)). 

Given a subspace A as above, we define the hyperbolic distance s{Pi, P2) between 
two points Pi — (zi, 2/1) and P2 — (z2, 2/2) in ^ Zi > 0, i = 1, 2 to be: 
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s{Pi,P2) := VTbi^7^-ln^2F+^^i-2^: if < zi, Z2 < 1 

otherwise it is defined to be equivalent to tlie standard euclidean metric. 

We define the parabolic hyperbolic distance between two points Pi = {zi,yi,ti) 
and P2 = {z2, j/2j ^2) with zi > 0, i = 1, 2 to be: 

s(A, A) := s(Pi, P2) + ^J\ti-t2\ 

where Pi = (zi,yi), P2 = (22,2/2)- 

Let < a < 1. We define the Holder space C"'^{A) in terms of the above 
distance. We start defining the Holder semi-norm: 

llfll - s„n \I{Pl)-I{P2)\ 

P^^P^eAn{{x..y)m'^:z>0] S[Pl,P2j 

and the norm 

ll/llcF(^):-|l/llco(^) + ||/lk.M) 

where ||/||co(^) := sup \f{P) \. 
PeA 

We say that a function f belongs to Cf'^'iA) if f° £ C°'{A°) and f G C^iA). The 

norm of / in the space C^'P{A) is defined as: 

\\f\\cf-''{A) '■= \\.f° WciA") + ll/llc'»(^)- 

Moreover, we define: 1 1 / | \co.p{a) — 1 1 /° I Ic«(^°) + 1 1 ./ I \c"(A) ■ 

Remark 1. We observe that /(w, y) e C"(^) if and only if /(z, y) e C"(^), where 
w ^ Inz. 

We say that a continuous function f on A belongs to C'^^f{A) if f° G C'^{A°) 
and f has continuous derivatives 

Jzj Jy} Jzz} Jzyj Jyy 

in the interior of A, such that 

2 \} ^ J )t ^ Jz,Z \ Jy ^ Jy ), Z Jzz, Z Jzy, Z [ Jyy — Jyy ) 

extend continuously up to the boundary. The norm of / in the space C'^^f{A) is 
defined as follows: 

2 2 

\\f\\c^^HA)--^\\J2^7f°\y'(Ao)+ E WDTD^fWcoa) 

m—O rn+n— 

Given f e C^+PiA), we say that f belongs to Cg+"'P(^) if 

r eC2 + "(^°) andzfzjy, Z^fzZ,ZfzyJyy 

extend continuously up to the boundary, and the extensions are Holder continuous 
on A of class C^'P{A). The norm of f in the space Cg '^'^{A) is defined as: 

2 

ii/iic.?+-(^):=iiriic^+"(.i=)+ E ii^'"^r^^/iic?-(-4) 

rn+n— 
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Remark 2. It follows by definition that fw — —pf + z^Pfz and fww — ~P fz + 
(1 — p) z^^'P fzz, which implies that: 

2 

r?2+n— 

Remark 3. The function / G Cg+"'P(^) if and only if /° G C2+"(^°) and / £ 
C^^"(^), therefore, the following norms are equivalent: 

l/llc,"+°-''(^) - ll/°llc2+°(^°) + ll/llc2+»(^)- 

Let T > 0. The definitions above can be naturally extended on the space-time 
domain Qt by using the parabolic distance ds^ = ds^ + \dt\. We define the space 
CfiQr) to be the standard Holder space with respect to the metric ds^ . We 
say that a continuous function f on Qt belongs to C'^'^^{Qt) if f has continuous 
derivatives 



JtT Jzt Jy: Jzzj fzy, fyy 

in the interior of Qt and f° has continuous derivatives that extend continuously 
up to the boundary and 

Z-P{f-f°), Z-P{ft-fn, z'-Pfz,Z~rfy, Z^~Pfzz, Z'-Pfzy, Z-Pifyy-f^y) 

extend continuously up to the boundary. The norm of / in the space C'^^p{Qt) is 
defined as follows: 

2 

ll/llc^+-=ll.rilc^+ E \\DiD^D{f\\co 

l+m+2j=0 

The function f belongs to CI'^°'^^{Qt) if f e C^+p{Qt), 

f,ft,zfz,fy and z^fzz,zfzy, fyy belong to Cg'^'iQT)- 



Throughout the paper k will denote a positive integer. We can extend these 
definitions to spaces of higher order derivatives. We denote by C^'P{Qt) the space 
of all functions / whose fc-th order derivatives D], D^ D\ f , i + j + 21 = k in the 
interior of Qt and z' D* D^yD\ {f — f° ), i+j + 21 = k exist and belong to the space 
C'iQT). We define C°^^p{Qt) = Hfe C^'PiQT)- 

We denote by C^+'^'PiQT) the space of all functions / S C^'PiQT) such that 
z'DlDiyDlf, iov i + j + 21 = k belong to the space C^^p^Qt)- The space 
C^'^°''P{Qt) is equipped with the norm: 

i+j+2l<k 

Remark 4. A function / G C^+°'^p{Qt) iff /° e C''+'^{Q^) and / e C''+°'{Qt)- 
Moreover, 

11/ IIcJ + °'P(Qt) — II f° llc'' + °.['=/21+-/2(Q5,) + 11./ |lc"» + °,[fc/21+°/2(Qj,)- 
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In the next paragraph we denote by Sq the half space x > in M^, by 5 the 
space S — So ^ [0, oo), and by St the space S x [0,T], for T > 0. The operator 
Lk ■■ C^+^+^'P^Qt) ^ C^+'^'pIQt) is defined as: 

(3.1) Lk[f] ■■= ft - {z'^aiifzz + 2zai2fzy + 0,22 fyy + blZ fz + 62/y + cf) 

where the coefficients {fly }i J are uniformly eUiptic and {fly, 6j, c} C C^^^^'^{Qt), 
{a22, 62, c} C C^+^+'--p{Qt), 1,3 = 1,2. 

Theorem 1. (Existence and Uniqueness) Let L^ be defined as above. Assume that 
<j> G C^+"'P(iS) and /o € Cf^+^+"'P(5o), and that 4>, /q are compactly supported in 
S and Sq, respectively. Then, for any T > 0, the initial value problem : 

(o 9^ / Lk [/] = (j) in St 

^ ' I f{;0) - /o on So 

admits a unique solution f £ C^+^+"'P(iSt). Moreover 

(3-3) ||/||pfc+2+a,p(^^^^ < C{T) (^||/o|lcj+2+-.p(5„) + ll0llcJ+°'P(5)j 

for some constant C'(T), depending on a, p, k and T . 

Proof. To solve the above Cauchy problem is equivalent to solve the following 

Cauchy problems (|3.4p and (|3.5p . 

The problem (|3.4[) is obtained by evaluating (|3.2p at z = and the problem (|3.5|) 

is obtained by solving the corresponding one for /. 

.04. / {Lk)o[f°]^<l>° in^x[Q,T] 

^ ' I r(-,0) = (/o)° onR 

(0 r^ f Lk [/] =0 in St 

^ ' I /(•,0) = /o on^o 

where the operators {Lk)o and Lk are defined respectively as follows: 
{Lk)o{f°)^al^f°y + hlf; + c°f°- 



(3.6) La; [/] ^ ft- {ail f WW + 2 ai2fwy + 0-22 fyy + bl fw + b2 fy + c f + G) 

with 



aij{w,y,t) 

bi{w,y,t) 

b2{w,y,t) 



■= a,j{x,y,t) 

:= {2p-l)aii{w,y,t) + bi{x,y,t) 

:= b2{x,y,t) 

c{w,y,t) := e^P''[p'^aii{x,y,t) -2pai2{w,y,t)+pbi{x,y,t)] 
G{w,y,t) := b2{w,y,t)gl{y,t)+a22{w,y,t)gly{y,t). 

By the assumptions on the operator Lk it is clear that the coefficients of the 
two operators {Lk)o and Lk satisfy classical conditions. We first find the solution 
f° to p.4|) . then we solve ()3.5p . By classical theory both problems have a unique 
solution. The function / defined by f{w, y, t) :— f°{y, t) + zP f{w, y, t) is a solution 
to (13:211. 
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Let (7'^+" and C'^+'^+°' denote classical parabolic Holder spaces. Then the fol- 
lowing inequalities hold: 

|/°||cfc+2+=(E+x[0,T]) < C'(r) (ll/o l|c'= + 2+a(B+) + Il5°llc'' + °(R+)) 
||/|lc;fc+2+c(5y) < C(T)[\\fo\\^k+2+c(So) + \\<P\\c'' + ''(S)) 

It follows that the solution to p.2p is unique and satisfies the inequality (|3.3p . D 

Next, we define the boxes in which we prove the Schauder estimates. Let < r < 1. 
We denote by Br{P) the box 

{/z\ z > 0, |a; — xq\ < e^ 
\V\-- |y-2/o|<r 
\t) to-r^ <t<to 

around the point P = \ Vo \ ■ We set Br to be the box around the point P = 
\toJ \h 

Remark 5. The choice of the box Br is made so that it has the right rescaling. The 
operators Lko and Lk are well understood on the corresponding boxes B" and Br- 

Theorem 2. ( Schauder Estimate) Assume that all the coefficients of the operator 

L f = ft - {z^ ailfzz + 2 Z ai2fzy + a22fyy + Z hfx + b2fy + cf) 

belong to the space C^^°'{Bi) and that the coefficients a22,62 o,nd c belong to 
(^fc+Q,p^Q^'j jpj. some numbers a, p in Q<p<l,0<a<l and satisfy 

«»jff >A|ep, vee7^n{o} A>0 

with ||ay|lcfe+°(Qy): Il^i|lcj+°(Q^)> ll«22|lcJ+°.P(Qy)JI^2|lcJ + °.!'(Q^)-l|c|lcJ+°.P(QT) - ^^ ' 

Then, there exists a constant C depending only on a, A and p such that 

for all functions f e C^+"-P{Bi). 

Proof. The proof follows by the same argument as in Theorem [1] and classical 
Schauder estimates for strictly parabolic operators. D 

4. The Degenerate Equation on the disc 

We will extend in this section the existence and uniqueness the Theorem [T] to 
the following class of linear degenerate equations: 

Lw := Wt — ( a^^Wij + b'' Wi + cw) 

on the cylinder V x [0,T), T > 0, where V denotes the unit disk in R^. The 
sub-indices i,j G {x,y} denote differentiation with respect to the space variables 
x,y. The matrix {a''-'} is assumed to be symmetric. Certain assumptions on the 
coefficients will be made so that this class of equations includes, under appropriate 
change of coordinates, the operator L given by (|3.ip . 
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We define the distance function s in 2? as follows: in the interior of V, s it is 
equivalent to the standard euclidean distance, while around any boundary point 
P G dV, s is defined as the pull back of the distance function induced by the metric 

T— 2 7 2 

as = -^r- + ay 

on the half space Sq = {{z, y) : z > 0}, via a map (p : Sq O V ^ V that flattens the 
boundary of the disk V near P. 

The parabolic distance is deflned by 



u r \t 



s(Pi, P2) + Vl*i-^2|, P1.P2 ev, < ti < t2 



We define the spaces C^+"'P{V) and C^+^+°''P{V). For a fixed smaU number S in 
< (5 < 1, we write 

v = v,_s/2 i){\JiVs{Pi)nv)) 
I 

for finite many points Pi € dV, I e /, with T>i_s/2 denoting the disk centered at 
the origin of radius 1 — 5/2 and 'Ds{Pi) denoting the disk of radius 5 centered at P;. 
We denote by V-i- the half disk 

V+^{{x,y)eV: x>0}. 

We can choose charts T/ : P+ -^ T>s{xi) DV which flatten the boundary of "D and 
such that T/(0) = Pi, I E I. Let {ip, ipi} be a partition of unity subordinated to 
the cover 

{v,_s/2AT^s{Pi)nv)} 

of P, with I g /. 

We define C^+"'P(X') to be the space of all functions w on V such that w £ 
C'=+"(2?i_5/2) and woTie C^+°''P{V+) for aU I e I. 

Also, we define C^+'^+"'P{T>) to be the space of all functions w on V such that 
w e C'=+2+°(2?i_5/2) and w o iPi g C^+^+"^p{V+) for aU I e /. Here C''+°' and 
(jk+2+a denote the regular Holder Spaces, while C^+°''P{V+) and C^+^+'^'P{V+) 
denote the Holder Spaces defined in section [3l One can show that both spaces 
Ck+a'P{V) and C^+'^+°''P{V) are Banach Spaces under the norms 



; 



and 



lkllcJ + 2+ = .P(I5) = \Hw\\c''+^+-{V,_s/2) +Y1 \\i^liwoTi)\\^k + 2+^,pf^^y 



I 
The above definitions can be extended in a straight forward manner to the par- 
abolic spaces C^^PiQ) and C2+"'P(Q) where Q is the cylinder Q = V x [0,T], for 
some T > 0. Before we state the main result in this section, we will give the 
assumptions on the coefficients of the equation 

Wt — a^-' Wij +b^ Wi + cw 

on the cylinder Q ^V x [0,T), i,j = 1,2. 

We first assume that for any (5 in < (5 < 1, the coefficients {a^-'}, {&*} and c 
belong to the Holder class C"(2?i_5/2 x [0, T]), which means that the coefficients 
are of the class C" in the interior of V. For a number 5 in < (5 < 1, let T; : 
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T>+ -^ VsiPi) n T) be the collection of charts which flatten the boundary of P, 
considered above. We assume that there exists a number S so that for every / S /, 
the coordinate change introduced by each of the T/ transforms the operator 

(4.1) L[w] = Wt — { a'-' Wij +b'^ Wi + cw ) 
on VsiPi) n P, into an operator Li on I?+ of the form 

Li [w] ^ wt - {x^ ail Wxx + 2xai2 w^y + 0-22 Wyy + xbiWx + b2Wy + cw) 

with the coefficients 5,^-, hi and c belonging to the class C^'^"(I?+), with 022, &2 
and c e C^+'^-P such that: 

for some number A > 0. 

We need the next Lemma to prove the invertibility of the operator L: 

Lemma 1. (Holder Interpolation) . For every e > there exists a constant C(e) 
depending on e, p, k and a such that for any g G C^^'^^°''P{Qs), the following 
inequality holds: 

(4.2) W^Dg I |pj+o.P(Q^) < e 1 1 .9 I IcJ+'+^'^CQ,) + C'(^) 1 1 9 I lc"=.p(Q,) • 
where "i? behaves like distance to the boundary. 

Proof. It follows by standard arguments. D 

The following existence result readily follows from Theorem [1] and the above 
discussion. 

Theorem 3. Assume that the operator L satisfies all the above conditions on the 
cylinder Q = T) x [Q,T]. Then, given any function wq G C^+^+"'^(I?) and any 
function g G C^ "'^(Q) there exists a unique solution w G Cs^^^"'^(Qt) of the 
initial value problem 

Lw — g in Q 

w(-, 0) — Wq on V 
satisfying 

(4.3) ||w||pfc+2+c,p(Q^ < C{T) (^||wo||cJ+2+°.pcD) + ll5llcJ+".p(Q)j 

The constant C{T) depends only on the numbers a, k, A and T . 

Proof. We can assume, without loss of generality, that wq = Q and that g is a 
function in C^'^'^'P {Qt) , which vanishes at i = 0. 

For 5>Q,&eiQ5 = Vx [0,(5] and denote by C';'^o^+°'-'P{Qs) and C,*^,^"'^(Q5) the 
subspaces of C^ +^+" {Qs) and C^ +" [Qs) respectively, consisting out of all functions 
which vanish identically at t = 0. Also, we denote by / the identity operator on 
CgQ''^{Qs). We will show that, if 5 is sufficiently small, there exists an operator 
M : C^+"'P(Q,) ^ C,"^+'+"'P(Q5) such that 

II LM-/ II < -. 
II 11-2 



This impfies that the operator L M : C^ a°''^{Qs) — > Cs o^'^iQs) is invertible and 

■■s,0 \S!&) ~* ^s,0 



therefore L : Cjg ^"'^{Qs) — > C^^'^'^iQs) is onto, as desired. 
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We begin by expressing the compact domain 2? as the finite union 

i>i 
of compact domains in such a way that 

dist(Po,9P) > ^ >0 

and for all/ > 1 

Vi = Bpixi)nv 

with Bp{xi) denoting the ball centered at xi E dV of radius p > 0. The number 
p > will be determined later. 

The operator L is non-degenerate when restricted on the interior domain I?o- 
Therefore, the classical Schauder theory for linear parabolic equations implies that 
L is invertible when restricted on functions which vanish outside 2?o- 

We denote by Mq : C,*^+"'^(r>o x [0,(5]) -^ C^^^+"'P{Vo x [0,(5]) the inverse of 
the operator L restricted on Vg. Next, we consider the domains Vi, I > 1, close to 
the boundary of V, which can be chosen in such a way that the sets Bp/4{xi) n P 
are disjoint. Denoting by B the half unit ball 

B = {(.T,y)eBi(0); x>0} 

and by Qg the cylinder 

Qs = Bx[0,d] 
we select smooth charts Ti : B ^ T>i, which flatten the boundary of V, i.e., they 
map Br\{x — 0} onto ViCidT) and have T;(0) = xi. This is possible if the number p 
is chosen sufficiently small. Under the change of coordinates induced by the charts 
Tj, the operator L, restricted on each Di x [0, S], is transformed to an operator L/ 
of the form 

Li [w] = wt- {x^ a]^wii +2x aj'^wi2 + apW22 + xbj wi +bf W2 + Q iu) 

defined on i? x [0, S]. Moreover, the charts T; can be chosen appropriately so that 
the coefficients of Li satisfy 

a^^^^J > A|^|2 > V^ G m2 \ {0} 

and 

for some positive constant A. 

Each of the operators L/ has the form of the model operators previously studied. 
Denote by 5*0 the half space a; > in R^ and by Ss the space Sq x [0,6]. Also, 

consider the subspace C^q ' (Ss) of C^~^°''^{Ss), consisting out of functions which 
are compactly supported on Ss- Then, Theorem [1] implies that for every / = 1, 2, ... 
there is an operator Mi : c'l'^o'^'" (Ss) -^ C^;^'^^°'^^{Ss) such that 

LiMi=I 

with / denoting the identity operator on C^ q ' (Ss)- Let Mi be the pull back of 
the operator Af/ via the chart T;. Next, choose a nonnegative partition of unity 
4>i; I — 0, 1, ... subordinated to the cover Vi; I — 0, 1, ... of 2? and also choose, for 
each / > 0, nonnegative, smooth bump functions ipi, < tjji < 1, supported in P; 
with ipi = I on the support of ipi. Then J2i>o 4'i = ^ ^-i^d V"; 4*1 = 4'i for all I. 
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We aim to show that the operator M : Cg^"'^{Qs) -^ Cj, o ^"'^(Qs) defined as 

satisfies 

\\LMg - 5||pj+..p(Q^^ < - ||g|lcj+=.P(Q^) Vg £ C^^a'^iQs) 
if the cover {P/} and S are chosen appropriately. Indeed, we can write 

LMg-g^J2^ ^i^'^i ^i9 " ^ 0/5 = $] ^i (LMi - I) ^ig + ^ [i, i/;, ] Mi (jug 
III I 

with [ L, ^i ] denoting the commutator of L and "0; . The commutator [ L, t/i; ] is only 
of first order and it can be estimated as 

Let e > 0. It follows via the Holder spaces interpolation from Lemma [T] that 

However, for each k we have 

\\Mi(t)i <?||pj+2+„,p(Q^^ < C||g||(^fc,=,p(Q^^ 
and therefore, since Mi(f>ig = at i = 0, 

|Af/0/5||pj,P(Q^) < C(5||5||p.+„,P(Q_^). 
It follows that if we choose S sufficiently small: 

On the other hand we have {LMo — I)(po9 — 0, while for / > 1, we can make the 
norm of each of the operators LMi — I arbitrarily close to zero by choosing the 
diameters of the domains Vi sufficiently small: 

II J2'^dLMi~I)(l3ig\\^,+..p^Q^^ < - ||5|lcj+=.p(Q_,) 

for all g S Cg^°''P{Qs), if p and S are both sufficiently small. The above inequalities 
give: 

WLMg-gWci+'^.P^Q^f < 2 ll5llcj+°.p(Q,) 

for all g G Csm^'^^Q^)- We conclude that for every g g Cg^""^{Qs) there exists a 
function w £ Csf^^^iQs) such that Lw = g. In addition 

(4.4) ||«;||c.j+2+o,p(Q^) < C||5||c.j+«.P(Q^) 

with C depending only on V and the constants a, fc, A and T. 

The last inequality implies we can extend the solution on a bigger interval. 
Hence, one can show that 

l^llcJ + 2+"-P(Q) - ^(^) (ll^o|lcJ+2+°.PCD) + llffllcJ + "-P(Q)j 

where the constant C(T) depends only on the numbers a, fc, A and T. 

This last inequality implies uniqueness. D 
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Finally, the following existence result follows from Theorem [3] and the Inverse 
Function Theorem between Banach spaces along the line of the proof of Theorem 
7.3 in [9J: 

Theorem 4. Let wq be a function in C^'+^+"'P(X'). Assume that the linearization 
DM{'w) of the fully-nonlinear operator 

Mw = wt — F{t, u, V, w, Dw, D w) 

defined on the cylinder Q = V x [0,T], satisfies the hypotheses of Theorem\^ at all 
points w G C^'+^+"'^(Q), with \\w — u'o|Ufc+2+c«,p,QS < /i, for some /i > 0. Th 



en, 



there exists a number tq in < to < T depending on the constants a, p, k, X and 
fj,, for which the initial value problem 

Wt — F{t,u,v,w, Dw, D^w) in 2? x [0, rg] 

w{-,0) —Wo on V 

admits a solution w in the space (j^+'^+<^,p(x> x [0, tq]). Moreover, 

for some positive constant C which depends only on a, p, k, A and /i. 

5. Global change of coordinates and existence in (7^+2+",? 

In this section we introduce a global change of coordinates which transforms the 
HMCF for a surface Eq into a fully-nonlinear degenerate parabolic PDE on V. 

Let 5* be a smooth surface close to S2- Let 5 : 2? ^ R'^, indicate a parameteri- 
zation of S on the unit disk T). 

S{u,v) — {x,y,z) &R^ which maps dV onto 5 Pi {x = 0} 

We denote by x„, x^,, Xum Xuv, x^ the partial derivatives of x with respect to u and 
V. The same notation will be used for the function y. 

Let r/ > be sufficiently small. Let T — {Ti,T2,T3) be a smooth vector field 
transverse to S. We define the global change of coordinates ^ : V x [—rj, 77] -^ M.^ 

by 



(5.1) \y] ^<^\ V \ =Siu,v) + wT{u,v) 
or more explicitly 

(5.2) X = Si{u,v) + wTi{u,v) 

(5.3) y ^ S2{u,v) +wT2{u,v) 

(5.4) z ~ Sz{u,v) A-wT'ilu^v) 
Let (5 > be small number, such that 

T3(u,w) = on V\Vi^s 

denoting, as above, the transverse vector field to the surface S. Notice that 
by choosing the smooth surface sufficiently close to the surface z = h{x,y), we 
can make S to depend only on the constant A which depends on the initial non- 
degeneracy conditions on the surface Sq. 
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We write the first and second derivatives of z with respect to x, y and t in terms 
of the first and second derivatives of w with respect to u, v and t. 
\i z = ho{x,y) then we compute the first and second partial derivatives of z with 
respect to x and y in terms of w = l{u, v) seen as functions of u and v. 

Let A be the Jacobian matrix relative to the transformation of coordinates: 

/ du dv\ 

dx dx 
du dv 

\dy dyj 

Let Vz, Vu, Vw be, respectively the gradients of z, u and v. 

-^(|)^(s)v..(|).(;)v.,.;r 

We denote by D^u and D^v the following matrices: 



A 



a b 
c d 



( 



D^u 



-\ 



dx"^ dxdy f '^i ci 

d'^u d'^u V "^i ^^2 

V dxdy dy"^ J 

/ d'^v d'^v \ 

d^ dxdy _ fh di 
d^v d^v ^ \di d2 
\ dxdy dy"^ I 
Based on the above, we define a2 := C\ , hi := di and denote ei 
to be the basis vectors. Let A^^ be the inverse matrix of A: 



D^v 



(1,0) 62 = (0,1) 



A-':= 



Xu, Xy\ / Xl X2^ 

Vu yvj \yi y2, 

Next, we introduce the matrices Bi and B2 which denote, respectively, the deriv- 
ative of the inverse matrix of A and A^^ with respect to x and y; Bi — — — 

B2 



dx ' 



dA 



Q — which can be computed as: 



axii 
axi2 



-bxi2 ayii+byi2 

-bx22 a yi2 + b j/22 

' CX11+ dxi2 cyii+dyi2 

CXi2+dx22 Cyi2+dy22 

Note that we are using the following notation: 



B2 = 



a;ii = Xy,u, a:;i2 = Xuv, 2:22 = Xyy, yn = y„„, yi2 = y™, 2/22 = yw 
The coefficients of the matrices D^u and D^v are evaluated as follows: 



= -A-Bi- Vu; 



= -A- Bi-Vv 



= -A- Bo- Vu; 



= -A- B2-Vv 
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dw 

Since: Vz — A ■ [ " ] + Zy^ [ §^ ] , wc obtain that: 

. Z^) 



ox \ZvJ \azuv+ozvv + Zy^-Q^ J \dyJ \'dl(hj, 



dy 



Vz = Ay 



-A- 



^ '^uu I ^ ^uv ' Zy 

C Zf^y -j- d Zyy -{- Zy 



dw 



dw \^\CZuw-\-dZyw) I dw I +^M) I g2 



dy 



d'w 



ff^w 
Wy 



The gradient of the function w, as well as its partial derivatives, can be expressed 
by using the matrix A: 



j§\ _ I awu + bwv 

CWu + dwy 



^1 - C aiwu + a{awuu + bwuv) + biWy +b(awuy + bwyy) 
■§^1 \ciWu + c{awuu +bwuy) + diWy +d{awuv + bwyy) 



d^w 



d^&y \ _ f a2Wu + a{cwuu + dwuv) +b2Wy +b{cwuy + dwyy) 



■jj^ / \C2Wu + c{cWuu + dWuv) + d2Wy + d{cWuv + dWyy) 

Q2y / \ ^ 



Using the substitutions above we get: 



—5— = A\^ wii+A\2 W12+AI2 W22+-S1 W1+B2 'W2+BI2 wi W2+BI^ WJ+B22 wl+Ci 
O'^x 



—5- = All W11+A12 W12+A22 W22+B', W1+B2 W2 + B12 Wi W2+-B11 Wi +i?22 W2+C2 



d^z 
dxdy 



^11 ^11+71^2 '"^12+^22 W22+B1 ■W1+B2 W2 + B12 Wi ■W2+B11 ■W1+B22 W2+C^ 



where the coefficients A'j; ,, Bf, C* with fc = 0, 1, 2, i,j = 1, 2 are defined as follows: 



41 
^11 

^22 

^12 

42 

42 
^22 

42 
^12 



== -a2 ( -Zyj + by,iy{Zu + Wl Zyj) + dy,iy{Zy + W2Zyj)) 

= -^2 {-Zyj + byw{zu + wi Zw) + dyw{zv + ■W2Zw)) 

= -2ab{-Zw + byw{zu + wi Zw) + dyw{zv + W2Zw)) 

== c^{zw - dyw{zy + W2Zw)) 

= d^iZyj - dyniiZy +W2Zyj)) 

= -2cd{~Zw + dyw{zy + W2Zn,)) 
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^?1 
^12 

-°ii 
RO 

S?2 



B'2 



C'> 



Bh 
Bh 
Bh 
Bl, 
Bh 
BI2 



Bl 



Ci: = 



Bl 
Bl 



C2 



-ac{-Zw + byw{zu + wiz^) + dyw{zy + W2^u.)) 

-bd(-Zw + hy^izu + wiz^) + dyw{zy + W2Z^)) 

-b(2acyuw + bcy^w + adyy^^)z^ 

-b(2acyuw + bcy^^ + adyyyj)Zy, 

-dibcyuw + adyuw + 2bdyy^)z^ 

-{{b'^c + a{b + 2c)d)yuw + d{b{2b + c) + ad)yyw)zn 



a \CXiiii ~r ^^uv )Zw '^x'^^Uvw^u ~r CCtyy^Zy CZyiy -r OCyuyZyj^ 



dzv 



b{2cyu 



ZZyiuJ 



OCtyvv^w ) 0,\ ^CZiiw \ ^CCLt/uw^v \ Ci IJvw^v 
Cy'Xuv \ IJuujZw I ^y^uv t yuvj^w)) 

yCyiiw "1" '^^yvw j^u ' Oy CZuw 1 ^\^yuw^u 1 Cyuw^v \ ^Ciyvw ^v 
[C[C -\- CljX^y + Cdy^iy -)- u. yvv }^w ) ClyC X^^^Zyj -r dy ^uw ~r CyXuv ~r yuuj^w 1 0,[yuw ^v 1 yuv^w ) ) ) 

(1 yCXuu ^^ d-^uv )^u Q-l^Ol^Cl^X-fj,^ -J- yiitlj ~r U^X^-^; -p yuv ) )^u ^^uu ~r C X-j^^Z^n" 
CdyXii^y ~r yuuj^v \ dy 'Z^^ -p dy^yZyjj OyOyCy^y -\- Ciyvv )^u C^uv ~r ^ -^uu^tJ ~r C 

^V-^Ml! ~r yUvj^V \ dyClyyyZy ^^ ^VV ) ) 

-2ab{ayuw + by^w) 

-2bd{ayuw + by^y,) 

-2(6^ + ad){ayuw + ^yi;«)) z^, 



Ad [C yuw ~r Ct yvw } ^w 

zcdycyuw 'T~ dyyuj jZiu 

ZdydyO yuw ^u ^uw ~r dyuw^v } ~r OyOyyuj ^u ~r dy^^yj Zy ^vw } } 

{d^ Xyy + d'^b{2XuV -\-yUu) -\-db^{XuV + 2?/^,^;) + B^ yyy) Zyj 

d yCXyy -\- d yuuj^W ^dOyOyyyjZy ^UW \ dy^yjZy "T CXyyZyj ~\- dyyyZyj j 

yZOyyyjZy "T ZdyyyjZy ' Z Zyy} "T" CXyyZyj "T dyyyZy^j 

d XyyZy d [0[ZXyy "t" yuUj^U ~^ ^UU I C-^UU^V "T dyyU^V } 

dOyOyXyy -]- Zyyy jZy ~~ ZZyy "T Z^CXyy ~T~ dy^y J Zy J \OyyyZy "r CXyyZy "T dyyyZy Zyy j 

ZCyyCyZyyj — dy^yjZy j "T d\ dyyyjZy -|- Zyyj j j G(2 Xyy Zyj 

Zdy CZyyj -\- CdyyyjZy "T d yyW^V " dZyyj j yC X yy H" C 

dyZXyy -J- yuu) I Cd yXyy + Zyyy j "+- Ct yVVj^W 

(2(1 XyyZy C XyyZy "T C yZyy dyZXyy "(" yuuj^V )\ 

CdyZZyy dyXyy "T ^yUVj^v) I ft (^ dyyyZy "T ^^^^ j 
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Evolution ofwt- The evolution equation of w is : 

(5.5) wt = zt 

where the function z satisfies the non Unear PDE: 

(5.6) Zt ^""^^^ " ^'^ 



I, J- I ZyjZxx ^ ^x^yZ^xy \^ ' ^x)^yy 

By replacing the partial derivatives of z in terms of the derivatives of w, we find 
the linearization L to the fully non linear operator given by the Equation 15.51 
Linearization. Let lu be a point close to the initial data Wq- Then, the linearization 
of the operator given by the Equation 15.51 is given by: 

L{w) = an wii + 2 ai2 wu + 022 W22 + biwi + b2'W2 + cw + d 

where its coefficients behave as follows: 

2 
an «;^ 511; ai2~^^gi2 

an « ^522; 61 « S^iJ^i^ /ii 

b2^^h2 

where {gij}, hi, c and d with i,j = f,2 are functions which belong to the space 
C". In particular 322 and 62 belong to Cf 'P. It follows that the coefficients of the 
operator L satisfy the same conditions of the operator of Theorem [S] 
We say that Sq ^s of Holder class (7^+2+". p if the function wq obtained belongs the 
space C^+^+°'^P{V). 

Based on the above definition, we state the following theorem: 

Theorem 5. Assume that the initial surface Eg belongs to the class (7^+2+".p ^^g^ 
satisfies the non- degeneracy condition (-k) . Then, under the coordinate change $, 
the HMCF with initial data the surface So converts into the initial value problem 

Mw = onVx[0, T] 
w — Wq at i = 

with Wo e C,*''+2+"'P(X') and 

Mw — Wt — F{t, u, V, w, Dw, D w) 

satisfying the hypotheses of Theorem^ 

As an immediate Corollary of Theorem |4] and Theorem [5] we obtain the following 
existence result: 

Theorem 5.1. Under the same assumptions as in Theorem [SJ there exists a number 
Tk > for which the HMCF with initial data the surface Sq admits a solution 
Ej on < t < r^,. Moreover, under the coordinate change $ the strictly convex 
part T,2{t) of Sj is converted to a function w{t) which belongs to the Holder class 
C^2+-.P(Q^),onQfe = Px[0,Tfe]. 

Theorem 6. Assume that the initial surface Sq satisfies the assumptions of Theo- 
rem]^ Then, the solution St of the HMCF is converted, via the coordinate change 
studied in Section\^ to a function w which belongs, for any positive integer k, to the 
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Hdlder class C^+^+°''P{Q), onQ ^Vx {0,T]. Moreover, for any t in < t < T 
we have 

(5.7) \\w\\^k+2+<..p^^^^^,j,^^ < Ca;(t, 1 1 Wo II cJ+° ■"(©)) 

Proof. We omit the proof as it is similar to the one done in [9 . Also, for more 
details we invite the reader to look at the Ph.D. thesis [5j of the first author. D 

6. The proof of the Main Theorem 

In this section we will give the proof of the Main Theorem stated in the Intro- 
duction. We will actually prove the following stronger result, where we relax the 
regularity assumptions on the initial surface. 

Theorem 7. Assume that the strictly convex part S2 of the initial surface Sq 

belongs to the class C^^"'^ and satisfies the non-degeneracy conditions (-k). Then, 

the HMCF 

dP K ^ , , 

with initial data the surface Sq admits a solution Ef which is smooth up to the 
interface, for < t < T. In particular, the interface Tt is a smooth curve for every 
< t < T which moves by the curve shortening Flow. 

Remark. It can be easily checked that if the initial surface satisfies the conditions 
of the Main Theorem, then it will satisfy the weaker conditions of Theorem [T] 

Proof. Assume that the strictly convex part S2 of the initial surface Sq belongs 
to the class C^+^'P and satisfies the non-degeneracy condition (*), then wc have 
proven existence for the HMCF in Theorem [51 

From Theorem[5]we have that w S Cg"'"^+"'^(Px (0, T]), for all nonnegative integers 
k. In particular this implies that for all integers k we have w{t) G C'^'^(P, for all r 
in < t < T. It follows that w is C°°'P smooth up to the boundary of V. Going 
back to the original coordinates, we conclude that the strictly convex part of the 
surface Et, < t < T is smooth up to z = and that the interface Ft is smooth. 

D 

7. Comparison principle 

In this final section we will give the proof of the comparison principle for the 
HMCF and we will show that the solution given in the Main Theorem is a viscosity 
solution. 

Proposition 7.1. (Comparison principle) Let Sq be a surface of class C^^"'^ that 
satisfies condition (•), and let S"*" be a smooth, strictly convex surface containing 
So at time i = 0, then the surface Sj obtained by evolving Sq by the HMCF, is 
contained in the surface S^ obtained by evolving T,q by the HMCF up to the time 
of existence of Et. Analogously, if Eq contains a smooth, strictly convex surface 
E~ at time t = 0, then the surface Et contains the surface E^ obtained by evolving 
Eq by the HMCF up to the time of existence of Ef . 

Proof. We observe that by the the classical maximum principle the surfaces Et and 
E( cannot touch were they are both strictly convex. Next, we suppose that there 
exists a time i where they first touch at a point P, then this cannot happen in 
the interior of the flat side. Hence P belongs to the boundary of the flat side. 
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Suppose St has the flat side on the a; = plane, then the tangent plane to the 
surface at the point P would be contained in the x — plane. Now if the two 
regions touch, then, because they are of class C^, the tangent to E^ at P would 
be contained in the x — plane. But S^ is strictly convex, hence this would imply 
that a part of E^ is inside Sf which leads to a contradiction. The second part of 
the proof is straightforward since if S contains a smooth, strictly convex surface 
E~ at time t — 0, then the two surfaces cannot touch at the flat side of Sj because 
the flat side does not move in its normal direction. Once again, by the classical 
maximum principle for parabolic equations the two surfaces cannot touch where 
they are strictly convex either. D 

Corollary 1. (Viscosity solutions) Let h ^ N, < j < I. Let Eq be a convex 
surface of class C^'^ that satisfies the hypothesis of the Main Theorem. Then the 
solution Et given by Theorem^is a viscosity solution of class C'''^. Moreover, this 
solution is unique upon satisfying the non- degeneracy condition (•). 



Proof. It is a consequence of the comparison principle. 



D 
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